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In a recent Comment larXiv: 1108.3 399 , Dzuba and Flambaum have referred to the disagreement 
of the results of our latest calculations of Tl electron and scalar-pseudoscalar (S-PS) electric dipole 
moments (EDMs) and Cs parity non-conservation (PNC) with some other calculations. We have 
responded to all their points and also discussed the larger issues related to them. We have attempted 
to find the reasons for the disagreement between the results of our calculations and those of others. In 
particular, we have found that the two important reasons for the discrepancies between the Tl EDM 
calculations of Dzuba and Flambaum and ours are the different choice of single particle orbitals 
and the treatment of the valence-core correlation effects. We have demonstrated by numerical 
calculations that the V^^'^ orbitals used by Dzuba and Flambaum overestimate the Tl electron and 
S-PS EDMs at the Dirac-Fock level. The failure of their suggested consistency test as interpreted by 
us is explained for systems with strong correlation like Tl. Also, the importance of understanding 
the physics underlying different theories on which atomic EDM and PNC calculations are based and 
comparisons between them are emphasized. 



In their Comment [l| , Dzuba and Flambaum begin by 
pointing out that our result for the electron EDM en- 
hancement factor for Tl is in disagreement with their 
calculation Q as well as Liu and Kelly's [1]. They then 
go on to state with reference to our paper where this re- 
sult was reported that "This is one more paper by this 
group presenting results which disagree with all other 
calculations. Others include, e.g. the calculations of the 
Tl EDM due to scalar-pseudoscalar (SPS) CP-odd in- 
teraction and parity nonconservation (PNC) in Cs." In 
this Comment, they do not, however, make any men- 
tion of our Cs EDM calculations based on the relativistic 
coupled-cluster (RCC) theory [H which agrees well with 
their calculations Q. With reference to this agreement, 
Dzuba and Flambaum had stated in their paper 0] that 
"Our calculations are in good agreement with the most 
recent and accurate calculations". The reference that 
they had given for the most recent and accurate calcula- 
tions was our Cs EDM work Q . Furthermore, they do 
not refer to our Ba"*" and Ra"*" PNC papers @, 0] in their 
Comment. In a recent paper Dzuba and Flamlaaum Q 
have reported that our Ba+ and Ra"'' PNC results arc in 
agreement with those of their latest calculations. 

We feel it is necessary to put the above remarks of 
Dzuba and Flambaum on our Tl EDM and Cs PNC cal- 
culations in perspective. A large number of calculations 
of the electron EDM enhancement factor (EDM-EF) for 
Tl have been performed using a variety of methods over 
the past four decades. We have mentioned in our paper 
that the results of these calculations vary from —179 
to —1041. This is also evident from Dzuba and Flam- 
baum's paper 0]. The reason for this variation is the 
peculiarly large electron correlation effects in Tl EDM- 
EF. While it is true that our result does not agree with 
those of other calculations which were carried out us- 
ing approximations substantially different from ours, it 



is equally true that the result of Dzuba and Flambaum 
agrees well with only the calculation of Liu and Kelly 
[j] . As we have explained in our paper this agreement 
is fortuitous. In their paper 0], Dzuba and Flambaum 
simply state this agreement, without giving any reasons 
for it. In contrast, we have endeavoured to explain 
the discrepancies between our results and those of Liu 
and Kelly Q and Dzuba and Flambaum [3[. For any 
scientific problem, the mere agreement or disagreement 
between two results is far less important than the rea- 
sons underlying them. It is therefore imperative to try 
to understand these reasons when one is comparing the 
results of two different theoretical methods. As for the 
differences between Liu and Kelly's work and ours, we 
reiterate the following points that we had stated in our 
paper "The following approximations were made in 
the former work (Liu and Kelly) (i) An approximation 
to only the one electron part of the EDM Hamiltonian 
is considered, thereby neglecting the important contribu- 
tions partly from the DF potential and largely from the 
two electron Coulomb interaction, (ii) Only the linear 
terms and a few selected nonlinear terms have been used 
in the calculations. The CC equations have not been fully 
solved even at the CCSD level as a coupled-electron pair 
approximation has been used to solve for the quadratic 
terms that have been taken in the unperturbed doubles 
equation, (iii) A few selected triple excitations are in- 
cluded only in the unperturbed singles amplitude equa- 
tions, where as, several dominant triples terms are com- 
pletely ignored in the unperturbed doubles equations. 
Thus, the contribution of triples is taken into account 
in a nonstandard way. (iv) The inner core is frozen up 
to the 4s orbital for the calculation of the unperturbed 
amplitudes, where as, for solving the perturbed doubles 
equations, 4s, 4p and 4d orbitals are further frozen. Such 
an inconsistent treatment introduces uncontrollable er- 
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rors. In contrast to the above drawbacks, we consider all 
the nonlinear terms arising from the single and double 
excitations. In addition, we consider the leading triple 
excitations in both the unperturbed singles and doubles 
cluster equations. We solve the unperturbed and per- 
turbed, closed- and open-shell equations, self consistently 
in the framework of the relativistic CCSD(T) approach, 
taking into account the excitations from all the core elec- 
trons." 

It is evident from the above comparison that Dzuba 
and Flambaum's point in their Comment [H that Liu 
and Kelly have used the same relativistic coupled clus- 
ter method that we have used to calculate the EDM-EF 
for Tl is extremely specious. It makes very little sense 
for us to spend a few months to try and reproduce their 
result, particularly when we know that it is based on in- 
consistent approximations leading to cluster amplitudes 
that are partially unphysical. Given these shortcomings 
of Liu and Kelly's work, a better way of testing the reli- 
ability of our Tl EDM-EF calculation would be to com- 
pare the results of our calculated quantities like allowed 
electric dipole transition amplitudes and hyperfine con- 
stants that are related to the EDM enhancement factor 
with available experimental data. That isprecisely what 
we have done in Table IV of our paper The overall 
agreement of our calculations with measurements is bet- 
ter than that of Dzuba and Flambaum [|. Liu and Kelly 
[3| have not performed any of these calculations. We are 
surprised that Dzuba and Flambaum have not referred at 
all to the comparison of the calculations of the relevant 
Tl transition amplitude and hyperfine constants in their 
Comment While these properties for Tl are not as 
sensitive to electron correlation effects as the EDM-EF, 
it is necessary to calculate them in order to get a sense 
of their accuracies as they are related to the EDM-EF. 

With reference to our comments about the compari- 
son between Dzuba and Flambaum and our calculations 
in our paper , we would like to categorically state that 
they are not misleading. We had not made any claims 
about the P and T violating Hamiltonian that Dzuba and 
Flambaum had used. Rather, we had surmised the form 
of their Hamiltonian as they had not given any informa- 
tion about it in their paper. In order to clarify, we repro- 
duce the relevant excerpt from our paper: "It appears 
from the previous work of Dzuba and Flambaum that 
the P & T violating Hamiltonian used in [10], considers 
only the internal electric field due to the nucleus and not 
the electrons; i.e., the entire two body Coulomb potential 
is neglected." It is still not clear to us from Dzuba and 
Flambaum's Comment whether they have included 
the gradient of the exact or the approximate two-electron 
Coulomb potential-they have not given the explicit form 
of their Hamiltonian. 

Indeed, we did mention in our paper that the V^^^ 
core, virtual and valence orbitals of Dzuba and Flam- 
baum are highly contracted-plots of some of these or- 
bitals are given in Figs, (a)-(f). We should have worded 
it more carefully and emphasized that the V'^~^ valence 



(6P1/2) and the virtual (7s) orbitals are highly contracted 
at large r relative to their V^~^ counterparts as shown 
in Figs, (d) and (f). The EDM-EF for Tl is very sen- 
sitive to the choice of 6pi/2, 7s and 6s orbitals and it 
is clearly desirable to use V^~-^ orbitals for an accurate 
determination of this quantity as they exhibit the cor- 
rect physical behaviour, particularly at large r unlike the 
orbitals. 

Table I highlights the large discrepancies between the 
V^-^ and Dirac-Fock (DF) resuhs for Tl EDM- 

EF and related properties. In retrospect, it appears that 
this is an important reason for the disagreement between 
the results of Dzuba and Flambaum Q and our cal- 
culations. The YiCC resuh for Tl EDM-EF would 
be —560 if one assumes that it changes in the same pro- 
portion to the DF as it does in the V^~^ case. This 
of course is not a rigorous assumption. It would not be 
out of place to mention here that in view of our DF re- 
sults for the two kinds of orbitals, the agreement between 
the CI-HMBPT result of Dzuba and Flambaum Q with 
yN-s orbitals and that of the linearized RCC calcula- 
tion of Liu and Kelly Q with V^~^ orbitals does indeed 
appear to be fortuitous as we had remarked earlier based 
on other considerations. 



TABLE I: Various results for different orbitals obtained 
using V^~^ and V^"^ potentials in Tl. 



States/Transitions 



6pi/2 

6pi/2 



Enhancement factors due to de 
-422.02 -507.15 

S-PS ratio in unit given in [9] 
5.02 6M 

Energy difference in cm~^ 



GPl/2 
GPl/2 



GPl/2 

6pi/2 



6pi/2 
6pi/2 



6s 
7s 



7s 



6s 
7s 



107507.67 
22714.75 



60294.73 
67140.58 



Reduced El elements in au 
2.08 2.02 
2.05 1.26 

Matrix element of Hi^j-,.,, 



44.26 
10.37 



60.89 
29.61 



It is clear from Table I of our paper Q that there are 
delicate cancellations between various valence-core cor- 
relations that we have evaluated to all orders using the 
RCC method in the singles, doubles and partial triples 
approximation. Dzuba and Flambaum [3| have calcu- 
lated some of these effects (Brueckener pair correlation) 
mainly by second order many-body perturbation theory 
and have neglected some others like structural radiation 
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and different classes of higher order RCC terms. The 
combined contributions from the latter two effects to Tl 
EDM-EF are about 30; which is by no means negligible. 
The approximations used in Q to calculate the valence- 
core correlations will not be able to capture the cancel- 
lations in our all order RCC calculations. 

We now turn to the S-PS EDM calculations for Tl. 
Three relativistic many-body calculations have been car- 
ried out for this quantity [3, 0, Dzuba and Flam- 
baum's calculation a gree s with that of Martensson- 
Pendrill and Lindroth jlO| . even though the approxima- 
tions used in the calculations are quite different. The 
latter work considers one electron effects to all orders, 
but its treatment of the two electron correlation effects 
is rather approximate. It only takes into account the 
Brueckner pair correlation which is not calculated di- 
rectly, but rather it is estimated on the basis of its con- 
tribution to the electron EDM-EF of Tl from Hartley 
et al's work [ll|, which had been performed using the 
same method that Martensson-Pendrill and Lindroth had 
used in their S-PS EDM calculation Therefore, the 
latter calculation is somewhat less reliable than the for- 
mer. The S-PS EDM calculations of Dzuba and Flam- 
baum and Martensson-Pendrill and Lindroth agree, but 
the electron EDM-EF calculation of the former is more 
than three times larger than that of the latter. This is 
indeed mind boggling and therefore one cannot attach 
any importance to the agreement of the Tl S-PS EDM 
results of the two sets of authors. Dzuba and Flambaum 
have not given any explanation for this paradoxical situ- 
ation. We had employed the RCC method in the singles, 
doubles and partial triples approximation and it is 
therefore not at all surprising that our calculation does 
not agree with the other two calculations which are based 
on different approximations. 

Dzuba and Flambaum [l| begin the last paragraph of 
their Comment by asserting that our Tl EDM-EF 
and S-PS EDM Q calculations do not satisfy a sim- 
ple consistency test. However, they conclude by mak- 
ing a rather tentative statement that these two calcu- 
lations may be internally inconsistent. The consistency 
test in their own words [l| is as follows: "the ratio of 
the EDMs due to two operators must be approximately 
equal to the ratio of the sp single-electron matrix ele- 
ments of these operators. This is because only short dis- 
tances, where single-electron energies can be neglected, 
contribute to the single-electron matrix elements of the 
CP-odd operators." We do not fully understand what 
this test is supposed to mean. Which s-p single elec- 
tron matrix elements are they referring to? 6s-6pl/2 or 
7s-6pl/2 or something else? They have not been very 
specific. Perhaps the meaning would have been clearer if 
they had supplemented what they have written by math- 
ematical expressions. Towards the end of their Comment 
[H , Dzuba and Flambaum seem to suggest that our value 
for the ratio they have mentioned in their consistency test 
is 115 de/C^^ IOI8 ecm and not 89 de/C^^ lO^^ ecm, 
which they probably think is the correct value. The for- 



mer is the value of the ratio of the electron EDM-EF to 
the S-PS EDM for Tl that we have got at the RCC level 
0, and our value for the same ratio at the DF level 
is 84 de/C^^ 10^^ ecm. It therefore appears that the 
consistency test that they have referred to in the con- 
text of our calculations means our ratios of the values 
of the two EDMs for the RCC and DF cases should be 
approximately equal. In other words, the ratio of the 
of electron EDM-EFs of the RCC and DF calculations 
is approximately equal to the the ratios of the atomic 
EDM to the S-PS constant for the same two methods. 
If our interpretation is correct then it is straightforward 
to show that the test that they have proposed does not 
hold in general, and in particular, we shall explain why it 
fails for Tl. In order to do so, we begin with the unper- 
turbed and first order perturbed wave functions in the 
RCC approach which can be expressed as 0, Q 



'{1 + 5(°)}|$„) 



and (1) 
{tW(i + 40)) +41)}|<I>,) (2) 

respectively, and the electron EDM enhancement factor 
or in general the ratio of the atomic EDM to a CP vio- 
lating coupling constant is given by the exact expression 

7^ = |D|vI/W)H-(vE'W|D|vI/(o)) 

= 7^l + 7^2 + 7^3 + 714 + 715 + 7Zq, (3) 



t(o) 



gT< '^£)gr' > ^j-^j each term along with its 



with D 

hermitian conjugate {he) term is given with subscripts 
from 1 to 6 with their sequence. To define the single 
and double excitations, we use the subscripts 1 and 2 
for all the RCC operators (for the detail see 0,0). As 
mentioned in [3, [l^ , the core and virtual contributions 
at the DF level comes at the lowest order through the 

RCC terms ($„|:D*°^ T(i)|$^„) and {^^\D^^^ si^'^\<^^), re- 
spectively. Denoting 7li = 71^^ + 7^^°'■'■ and 7^2 = 
71^'' + 7lf'^'', the total DF resuh is given by 71 df = 
7Z^^ + 7Z^^ . With these notations, the above expres- 
sion can be rewritten as 



7Z 
7Zdf 



= 1 



7Zf 



7^3 



7^4 



7Zdf 7Zdf 7Zdf 7Zdf 
7I5 , T^6 



71 



DF 



71df' 



(4) 



All the terms on the right hand side of the above equa- 
tion are built out of different combinations of D, T^'^\ 
S'i"^ r(i\ and Si^\ The last two operators; i.e. T^^' 
and si^'^ contain one order of the EDM interaction and 
all orders in the residual Coulomb interaction. The elec- 
tron EDM interaction Hamiltonian is given by 



H 



EDM 



2icd,^ Pjj^p: 



5 2 

A^I^J IjPj 
3 



(5) 
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and the S-PS interaction Hamiltonian is given by 

H'edm = '-^CsY.Pn'.P'nuAr). (6) 

where de is the intrinsic e-EDM, 7^ is a pseudo-scalar 
Dirac matrix, C5 is the dimensionless S-PS constant and 
Prmc(^) is the jth electron density over the nucleus. The 
corresponding single particle electron EDM and the S-PS 
matrix elements are given by 



^-P,{r)Q,{r) 







dPi dQj dQ, dPj 
dr dr dr dr 



(7) 



and 



dr{P,{r)Q.j{r) 



+Q:{r)P,ir))pi^Ar), 



(8) 



respectively, where I and / are the orbital quantum num- 
bers for the large and small components of the Dirac 
orbitals. It is clear that the above two matrix elements 
are numerically quite different, even though both of them 
are large in the nuclear region. The amplitudes for T^^\ 
and si^^ are therefore different for the two EDMs; par- 
ticularly in Tl where the electron correlation effects are 
strong. This is because these amplitudes represent dif- 
ferent processes involving the EDM and the correlation 
effects. Therefore „^ would in general be different for 
the electron EDM and the S-PS interaction making the 
consistency test that Dzuba and Flaumbaum have sug- 
gested invalid for systems with large correlation effects; 
although it could hold in situations when the correla- 
tion effects are weak or moderately strong. This is re- 
flected in the results of our RCC calculations of the two 
EDMs for Tl-see the relevant tables The trends 

exhibited by some of the RCC terms in the two cases 
are different-the most striking being the relative contri- 
butions of DT^^^ and DSi^\ As a result, the relative 



contributions of — and 



are quite different. 



We would now like to clarify certain issues related to 
our parity nonconservation (PNC) calculation in Cs that 
the Dzuba and Flambaum [l| have commented on. It 
is certainly not true that our Cs PNC result does not 
agree with the results of other calculations of this quan- 
tity. In an earlier work, we had reported our Cs PNC re- 
sult as 0.902(4) x 10~^^zeao— Qvv/-^ from a preliminary 



calculation [1^ and more recently we obtained 0.8892 x 



lQ~^'^ieaQ—Qw /N (triple excitations involving the va- 
lence and two core electrons were incorporated in the lat- 
est result) [11] using the RCC method. The latter result 
is in agreement with another recent calculation by Por- 
sev et al [ll|, who have got 0.8891 x IO~"zeao-Qiv/iV 
in the framework of the Dirac-Coulomb approximation 
using unsealed wave functions by a sum-over-states rela- 
tivistic coupled-cluster (CC) method at the singles, dou- 
bles and leading order triple excitations which is at par 
with our approach. However, Porsev et al have mini- 
mized the uncertainty in the calculation by scaling the 
wave functions and they have added other corrections to 
get the final result. They have used their final result to 
probe new physics. In our work, we have only demon- 
strated the important role of correlation effects using an 
ab initio RCC method and we have not put any effort 
to minimise the uncertainty at the Dirac-Coulomb level 
or considered higher order relativistic and nuclear correc- 
tions to probe new physics using our results. Our result 
is not in serious conflict with any other calculation if it 
is viewed in the right context. The reported results of 
Dzuba and Flambaum are 0.9001 x 10~^^ieao— Qw/iV 
[13 and 0.9078 x IQ-^^ieao-Qw /N [l^ at the Dirac- 
Coulomb approximation. The same authors had also re- 
ported the s-d PNC amplitude calculations in Ba"*" and 
Ra+ [13 about a decade ago [T^]- They had obtained 
2.17 and 2.37 (in xlO'^^ieao-Qw /N) for Ba+ and 42.9 
and 45.9 (in xlO^^^ieao— Qw/A^) for Ra+ using a mixed 
state and the sum-over-statcs approaches, respectively. 
In their Comment, the same authors had expressed the 
view that it is necessary to obtain accurate results for 
these calculations in order to use them for inferring new 
physics [l|. But their dual results could not have been 
useful for this purpose. Accurate PNC results for these 
ions were only identified when we employed the RCC 
method and reported our results as 2.46(2) and 46.4(1.4) 
(in -xlQ-^^ieao-Qw/N) for Ba+ and Ra+ [3, re- 
spectively. We learnt from a recent work of Dzuba and 
Flambaum that using one of their advanced methods, 
they were able to improve their results to 0.29 and 3.4 
(in X lQ~^'^ieaf)—Qw) [3 compared with our results 0.304 
and 3.33 (in x lO^^^ieao— Qiy) in Ba"*" and Ra"*", respec- 
tively. We quote below from their paper [^i "There is 
also good agreement with Sahoo et al. for Ba"*" [21] and 
with Wansbeek et al. for Ra+ [17] 

In conclusion, we have responded to all the points 
raised by Dzuba and Flambaum in their Comment [l| 
on our Tl EDM and Cs PNC calculations. The reasons 
for the discrepancies between Liu and Kelly [3] and our 
calculations which were discussed in our paper [^ have 
been reiterated. The two main reasons for the disagree- 
ment between Dzuba and Flaumbaum and our Tl EDM 
calculations have been identified as the choice of orbitals 
Dzuba and Flambaum and F^"^ by us) and 
the treatment of valence-core correlation. 

It has been pointed out that the the agreement of the 
Tl electron EDM results of Liu and Kelly and Dzuba and 
Flambaum is fortuitous. We have provided strong evi- 
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dence to show that the agreement of the Tl S-PS EDM 
results of Dzuba and Flambaum Q and Martensson- 
Pendrill and Lindroth raises many more questions 
than it answers. We have argued that the consistency 
test that Dzuba and Flambaum have referred to in their 
Comment involving the electron and S-PS EDMs as in- 
terpreted by us is flawed. Our RCC results for quantities 
that are associated with Tl EDM are in better agreement 
with experimental data than those of Dzuba and Flam- 
baum, which these two authors do not mention in their 
Comment. Contrary to the comments made by Dzuba 
and Flambaum on our Cs PNC work, our result does 
agree with that of another calculation that was performed 
at the same level of approximation [3l . 

Dzuba and Flambaum refer to three calculations (Tl 
electron EDM, Tl S-PS EDM and Cs PNC) in their Com- 
ment for which their results are different from ours. How- 
ever, in the same Comment, they do not mention the re- 



suhs of our Cs electron EDM, Cs S-PS EDM, Ba+ PNC 
and Ra+ PNC calculations which are in agreement with 
theirs. When the result of a calculation of ours agrees 
with theirs as in the case of Cs EDM it is described 
as accurate by them Q , but on the other hand if it does 
not agree then they make adverse remarks against it as 
exemplified in their Comment [l[ . Dzuba and Flambaum 
have not made any attempt to explain why a particular 
result agrees or disagrees with their result P, they 
have merely stated whether it agrees with their result 
or not. Such an approach provides no insights into the 
reasons for the agreement or disagreement between the 
two methods that were used to obtain these results. An 
understanding of the physics embodied by the methods 
used in the calculations of atomic EDMs and PNC and 
comparisons between them are prerequisites for making 
further progress in this field. 
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Piol ol' iLiriit: components til' win. lor Sdj^ 




R (a.u.) 



(a) 5ci3/2 orbital wave function 



Piol ol' iLiriit: components til' win. I'tir Silj/i 




R (a.u.) 




(c) 6si/2 orbital wave function 

FIG. 1: Comparison of magnitudes of the 5d^/2, 5d^/2, 
and 6si/2 single particle wave functions using V^^^ 
and V^^^ potentials. 
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Piol ol' iLiriit: components of win. lor 6pi 




R (a.u.) 



(d) 6pi/2 orbital wave function 



Plol ol' iLiriit: components til' win. I'or 6pi 




R (a.u.) 



(e) Qp'i/2 orbital wave function 



Plol ol' lariie componenl.s of wfn. lor 7S|/i 




R (a.u.) 



(f) 7si/2 orbital wave function 

FIG. 2: Comparison of magnitudes of the 6pi/2, 6^3/2, 
and 7si/2 single particle wave functions using V^~^ 
and V^^^ potentials. 



